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ABSTRACT 


A multiple-realization particle trajectory scheme has been developed and applied 
to the numerical prediction of confined turbulent fluid-particle flows. The example 
flows investigated include the vertical pipe upflow experimental data of Tsuji et al. 
and the experimental data of Leavitt for a coaxial jet flow, comprising a par- 
ticle-laden central jet and a clean annular jet, into a large recirculation 
chamber. The results obtained from the numerical scheme agree well with the 
experimental data lending confidence to the modeling approach. The multiple- 
realization particle trajectory turbulent flow modeling scheme is believed to 
be a more elegant and accurate approach to the extension of single-particle 
hydrodynamics to dilute multi-particle systems than the more commonly employed 
two-fluid modeling approach. It is also better able to incorporate additional 
force terms such as lift, virtual mass and Bassett history terms directly into 
the particle equation of motion as appropriate. This makes it a suitable can- 
didate for particle migration studies and an extension to situations involving 
liquid particulate phases with possible propulsion applications, such as in 
spray combustion, follows naturally. 
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INTRODUCTION 


Turbulent fluid-particle flows are encountered in numerous technological 
applications such as fluidized- bed combustors and pulverized coal gasifiers and 
combustors as well as in atmospheric studies involving the dispersion of pollutants. 
The modeling of such turbulent flows involving the presence of a dispersed phase 
made up of small, light particles further complicates the already complex phenomena 
encountered in single phase turbulent flows. However, the need to optimize the design 
process in technological applications involving turbulent fluid-particle flows or 
enhance the prediction accuracy of atmospheric dispersion models makes it impossible 
to avoid the quest for a deeper understanding of the fundamental problems. Besides, 
the various interacting complex phenomena encountered in the modeling of this class 
of flows offer a very rich source of challenges to the fluid flow researcher. 


The propulsion systems for space transportation vehicles, in particular the 
liquid-fueled variety, will benefit directly from an improvement in the modeling of 
turbulent fluid-particle flows. This is because such an improvement will translate to a 
better understanding of the mixing and combustion phenomena in spray combustion 
processes. Turbulent fluid-particle flows involving solid particles are simpler to model 
than fluid-droplet or fluid-bubble flows due to the added degrees of freedom in the 
latter associated with the deformation of the discrete entities of such a dispersed fluid 
phase. A study of turbulent fluid-solid particle flows is thus useful in elimi- 
nating the effects of the breakup or coalescence of droplets and bubbles from 
other particle-turbulence interactions encountered in such flows. 

The two common approaches adopted in the literature for the modeling of 
two-phase flows are the homogeneous and the separated models. The former is 
applicable to situations in which the mean slip between the phases is small and the 
design parameters of interest are of the bulk variety such as the pressure drop or 
mass fluxes. In situations where more detailed information about intra- or inter- phase 
behavior is of interest, or where there is substantial segregation of the phases, the 
separated two-phase models are invariably preferred. For such flows, another major 
decision has to be made with regard to the scheme for the description of the 
dispersed phase - whether to adopt an Eulerian or a Lagrangian approach. Important 
considerations necessary for deciding which approach to adopt include the concentra- 
tion of the dispersed phase which influences the mean separation distance between 
particles. The relative magnitudes of this length scale as well as the particle size and 
the microscale of the underlying turbulence in the continuous phase help to determine 
whether the dispersed phase can be treated as a continuum and thus described using 
the Eulerian approach or whether a Lagrangian description of the dispersed phase will 
be more appropriate. 

In the following, we present a discussion of turbulent fluid-particle flow 
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modeling in which the continuous phase is described using the continuum Eulerian 
approach while a Lagrangian description is adopted for the dispersed phase. We shall 
restrict ourselves to confined flows and thus include a discussion of the treatment of 
solid boundaries using the Eulerian - Lagrangian scheme. 


2. PARTICLE TRAJECTORY SCHEMES 


In the Eulerian - Lagrangian modeling of two-phase flows, the continuous fluid 
phase is described using the standard single phase continuum equations. How- 
ever, the dispersed phase is modeled by computing for individual particles the 
trajectories and temperature histories where appropriate. The dispersed phase 
velocity and temperature fields are subsequently obtained from information 
obtained from the realization of a sufficiently large ensemble of particle tra- 
jectories. 

The use of a particle trajectory scheme in the modeling of turbulent fluid- 
particle flows represents only a subset in the field of computer simulation using 
particles as discussed by Hockney and Eastwood [1981], Other important applications 
particle schemes discussed by Hockney and Eastwood include the modeling of covalent 
and ionic liquids, stellar and galaxy clusters, plasma and semiconductor devices. 

In fluid dynamic applications, the Particle-In-Cell (PIC) method of Harlow 
[1964] and, later, the Particle-Source-In Cell (PSI-Cell) method of Crowe et al. 
[1977] have received considerable attention. In the present investigation, the PSI-Cell 
method has been adopted as the basis for the Eulerian - Lagrangian model developed. 

The usual starting point for the development of fluid-particle flow theory is the 
consideration of the motion of a single particle in an infinite fluid. The nature of 
such a single-particle flow has been investigated by numerous researchers including 
Bassett [1888], Boussinesq [1903], Oseen [1927], Tchen [1947], Corrsin and Lumley 
[1956] .Hjelmfelt and Mockros [1966] and Maxey and Riley [1983] and is relatively 
well understood for flows both within and outside of the Stokes flow regime. In the 
Eulerian treatment of the dispersed phase, the single particle flow theory is adopted 
directly to describe a multi-particle system and the validity of such a step is assumed. 
However, in the Lagrangian particle tracking approach, the focus remains on single 
particle hydrodynamics for obtaining an ensemble of statistical realizations, in this 
case the particle trajectories, which are then analyzed using the well established 
mathematical theory of statistics to extract the required phase information. 


In the presence of turbulence, particle trajectories are not deterministic due to 
an imposition over the mean velocity of a rapidly fluctuating random velocity 
component. This additional velocity component due to turbulence enhances the 
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dispersion of the particles, in aggregate, while the presence of the particles in the 
continuous phase, even in relatively small concentrations [Al-Taweel and Landau 
1977], does modify the underlying turbulence appreciably. This ’two-way cou- 
pling’ between the turbulence and the particulate phase exercises considerable 
influence over the evolution of such flows. These important effects will be 
considered later. 


3. GOVERNING EQUATIONS 


The field equations for the continuous phase in the Eulerian - Lagrangian 
scheme are the same as those for single phase flows except for the addition of an 
extra ’source’ term which accounts for the influence of the particulate phase on the 
continuous phase. The equations are written in a generalized form as 

~ 3 *\ ( Pu i*)°ax7 ( reff fir) +s + s p (3.1) 

where u; are the instantaneous velocity components, T^ff the effective exchange 
coefficients, S the usual single-phase source terms, S p the source terms due to the 
particulate phase and <$> any of the field variables such as velocity component, 
temperature for flows involving energy exchange, turbulence kinetic energy or its 
dissipation rate. 

The simplified form of the particle trajectory equation in which only the 
hydrodynamic drag term between the phases is retained [Adeniji-Fashola and Chen, 
1987] is 


dv, _ (u j - vj) 

dt “ x* 


(3.2) 


where, in general, the fluid and particle velocities, uj and vj respectively are made 
up of a mean and a fluctuating component and T * is a particle response time 
defined in terms of the particle relaxation time t* which is valid for particle motion 
within the Stokes regime. Thus we have 


Uj = Uj +Uj' 

(3.3) 
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where 
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An expression similar to equation (3.2) for the particle temperature history 
written for a particle thermal equilibration time t-j-H can also be written for 
flows involving energy transfer [Chen and Adeniji-Fashola, 1987], 


4. PARTICLE-TURBULENCE INTERACTION 


A very important aspect of the modeling of turbulent fluid-particle flows 
is the particle-turbulence interaction problem. Turbulence kinetic energy 
extracted from the mean flow kinetic energy of the continuous phase is partly 
dissipated by the smallest eddies and partly imparted to the particles thus 

enhancing the dispersion of the particulate phase. This ’two-way coupling’ 

referred to earlier - modulation of the kinetic energy of turbulence by the 
particles and enhanced dispersion of the particles by the turbulence will now 
be discussed in a little more detail. It is pertinent to point out at this 

point that the turbulent dispersion phenomenon is primarily responsible for the 

considerable enhancement in mixing observed for turbulent flows when compared 
with laminar flows. 


TURBULENT DISPERSION 

The turbulent dispersion phenomenon is very closely related to the inter- 
action between individual particles and turbulent eddies. A particle normally 
interacts with a series of eddies as it moves through the fluid. The particle 
trajectory scheme attempts to simulate this interaction by tracking each repre- 
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sentative computational particle through a succession of turbulent eddies con- 
tained within the domain of interest. Figure 1 is a schematic illustration of 

this interaction between particle and eddies and in relation to the computa- 
tional cells. As discussed by Gosman and Ioannides [1981], a particle interacts 
with a given eddy for a period of time which is the minimum between an esti- 
mated particle transit time within the eddy, t tr and an eddy lifetime, t e . The 

particle transit time is obtained as the solution of the linearized equation of 
motion of the particle while the Lagrangian time scale of the turbulent eddy is 

obtained from length and velocity scales of the turbulence which are extracted 

from a k- turbulence model. Thus, 


tint 35 Min [t e ,t tr ] 


(4.1) 


where 


ttr® -t*ln [1.0-1 e/tjuj-v.l] 

and 


(4.2) 


te= ' e /(2k/3) 1/2 


(4.3) 


The eddy length macroscale, l e is defined in terms of the kinetic energy 
of the turbulence, k and its dissipation rate, £ as 



C 2/4 3/p 
k d /£ 


(4.4) 


In a stochastic formulation of the particle trajectory scheme which is the 
case in the present study, the fluctuating component of the fluid velocity, u\ 
is obtained from a Gaussian distribution of values having a zero mean and a 
standard deviation, Cjj given by 

G ii • (2k/ 3)^2 (4.5) 
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The Gaussian distribution is, however, not expected to be appropriate, in gen- 
eral, for describing non- homogeneous, non-turbulent flows. 


TURBULENCE MODULATION 

The presence of particles, even in very small concentrations, has the 
effect of modulating the turbulence intensity, the direction of modulation 
being influenced by the mean particle size and the level of modulation by the 
particle loading. This turbulence modulation effect was observed experimentally 
by Moderrass et al. [1984] and Tsuji et al. [1984] and attempts to mathemati- 
cally characterize the phenomenon include those of Al-Taweel and Landau [1977] 
and Chen and Wood [1985]. The interphase interaction force terms between par- 
ticles and the continuous phase are reflected as extra dissipation terms in the 
modeled equations for k and B when the former are included in the derivation of 
the field equations for the latter. The earlier attempts to implement these 
turbulence modulation models have been mostly within a two-fluid formulation in 
which the two phases are described as two interpenetrating continua viewed from 
an Eulerian framework. Equations (4.6) and (4.7) from Chen and Wood [1986] show 
the extra dissipation terms due to the turbulence modulation effect of the par- 
ticles for such a two-fluid formulation: 
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The term TH1 in equation (4.6) is the turbulence modulation term due to 
the mean slip while the terms TH2 and TH3 are due to the particle slip velocity 
at the fluctuating level. The model is valid for the situation 


t e > t* > tj^ , where tj^ = (V/6 )*/ 2 (4.8) 


is the Kolmogorov time scale. The model described above has been incorporated 
into the particle trajectory scheme of the present study. 



5. NUMERICAL SCHEME 


The set of governing differential equations describing the evolution of 
confined turbulent fluid particle flows cannot, in general, be solved analyti- 
cally thus requiring the adoption of a numerical procedure. For the continuous 
phase, the governing Eulerian equation set is solved using the SIMPLE algorithm 
of Patankar and Spalding [1972] and Patankar [1980], The overall scheme adopted 
for the solution of the governing equations is similar to that suggested by 

Crowe et al. [1977] and illustrated in Figure 2. An alternative scheme more suited to 
time-dependent flows was later presented by Dukowicz [1980] and further developed 
by Cloutman et al. [1982] and Amsden et al. [1985]. 

First, the "clean" fluid flow field is obtained by solving 
the continuous phase governing equations. This is done using a staggered grid 
distribution in which velocity cells are centered about the edges of the scalar 

cells. Next, particle trajectories are computed for a predetermined number of 
representative particles such that a statistically stationary solution is 
obtained for the overall particle flow field. The particle trajectories, and 
temperature history where appropriate, are obtained by solving for the particle 
the non-linear ordinary differential equations of motion and the energy equa- 

tion subject to the currently existing continuous fluid flow and temperature 
fields. A fourth order Runge-Kutta algorithm is used for this purpose. During 

the calculation of a particle’s trajectory and temperature history, the sources 
of momentum, energy, kinetic energy of turbulence and its dissipation rate, all 
due to the particle motion, are accumulated for each computational cell trav- 
ersed. The form of the source terms have already been presented elsewhere [Ade- 
niji-Fashola and Chen, 1987] and so will not be repeated here. These source 
terms are then used in the next global iteration on the continuous phase field 
equations until convergence is attained. It was found that source term relax- 
ation was required to achieve stability of the global iteration scheme for some 
of the example flow problems studied. 


PARTICLE SOURCE FIELD CONTINUITY 

A necessary condition to obtain a globally converged solution is to ensure 

the continuity of the source fields as was also pointed out by Durst et al. 

[1984], In order to ensure compliance with this important requirement, it is 

necessary to ensure the computation of source terms for each cell traversed by 
each computational particle through a judicious choice of the particle integra- 
tion time step as well as have particles start from as many locations as is 
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practicable within the relevant portion of the inlet plane. In the present 
study, particles are uniformly distributed in physical space at the appro- 
priate portion of the inlet plane of the computational domain in contrast to 
the scheme of Durst et al. [1984], in which particles are introduced only at 
grid nodes. The smooth profiles they obtained are very likely to be a 
consequence of the deterministic nature of the particle trajectories used in 
their study. 


INTEGRATION TIME STEP 

The choice of appropriate time steps for the integration of the particle 
equations of motion is very vital to obtaining a globally converged solution 
and smooth averaged particle flow fields. For the complex confined turbulent 
fluid-particle flow problems in general, some of the relevant time scales 
include the Lagrangian or macro time scale (eddy lifetime) of the turbulence, 
t e ; the Kolmogorov or the micro (dissipation) time scale of the turbulence, tj^; 
the particle relaxation time, t*; the particle residence time within a computa- 
tional cell or the whole computational domain tj^. Also relevant to the stochastic 
determination of the particle turbulent intensity are the particle transit time within an 
eddy, t te and the particle eddy interaction time, t| nt . The integration time step is 
selected to ensure adequate resolution with regard to the trajectory and temperature 
evolution while ensuring computational efficiency by avoiding unnecessarily small time 
steps. 


In the present study, a variable integration time step scheme was devised. 
An upper bound on the time step through any computational cell was imposed 

based on an estimated particle residence time for that cell and with the par- 
ticle being constrained to undergo about four integration steps within the 

cell. Without this restriction, the possibility of a particle overshooting one 
or more cells, possibly due to a sudden reduction in cell dimensions in a 

non-uniform grid domain, exists. Such a situation will result in a failure to 
compute the relevant source term contributions for a cell that was actually 
traversed by the particle. The consequence will be a lack of smoothness in the 
particle source distribution and, possibly, divergence of the global iter- 
ations. 

Also, for the reason of ensuring a smooth evolution of the particle tra- 
jectory and temperature history, a further restriction on the integration time 

step, At < t* , is made. The particle-eddy interaction time is determined and 
controlled independently of the integration time step. 
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PARTICLE AVERAGED PARAMETERS FROM PARTICLE TRAJECTORY STATISTICS 


One of the problems associated with the use of the Lagrangian particle 
trajectory approach, highlighted by Smoot and Smith [1985], is the difficulty 
of extracting smooth mean particle flow and temperature fields from the statis- 
tics of trajectories and temperature histories obtained for representative com- 
putational particles. In the present study, the fluid properties utilized in 
the particle trajectory and temperature history calculations are the linearly 
interpolated values in which the four nearest neighbors regarding the par- 
ticle’s current location are used, resulting in second order accuracy [Sirig- 
niano, 1983]. The details of the extraction of particle mean flow and tempera- 
ture fields information from the particle trajectory and temperature history 
statistics are available in Adeniji-Fashola et al. [1988]. 


BOUNDARY CONDITIONS 

The definition of a fluid flow problem becomes unique through the speci- 
fication of the boundary conditions after the governing differential equations 
are outlined and an appropriate closure of these equations is effected. The 
example flow problems investigated in the present study include vertical pipe 
upflow and horizontal recirculation chamber flow. However, rather than define 
the boundary conditions specific to each flow problem separately, the more 
efficient approach of defining generic boundary condition types is adopted. It 

then becomes a straightforward exercise to construct the boundary conditions 
for these and other specific flow situations of interest. 

Inlet Plane: 

The specification of the inlet plane boundary conditions for fluid flow 
problems is very important, as was discussed by Sturgess et al. [1983] and 

Westphal and Johnston [1984], since this influences significantly the 
subsequent evolution of the flow, especially in the case of parabolic flows for 
which the inlet plane conditions constitute the initial conditions for the 
solution of the governing differential equations. 

In order to correctly simulate a given fluid flow experiment numerically, 
the ideal specifications for the inlet flow variables are the actually measured 
values. The complete set of measured inlet flow variables is, however, hardly 
ever available. In the absence of such detailed experimental information, uni- 
form profiles are commonly specified for the axial velocity and temperature 

profiles of the continuous phase flow at the inlet plane. The turbulent kinetic 
energy is usually assumed to be a percentage, between 3 and 20%, of the inlet 
flow mean kinetic energy. The kinetic energy dissipation rate at the inlet is 
then obtained as 
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( (jV4 k 3/2 ) / , d 


(5.1) 


£ = 


where lj, the dissipation length scale, is specified as a fraction of the 
characteristic length scale at the inlet. 

For the particle trajectory and temperature history calculations, the 
initial velocity and temperature slip values relevant to the particular flow 

problems are employed in setting the required inlet conditions. 

Exit Plane: 

At the exit plane, the usual boundary condition imposed for any flow var- 
iable, <|> , is / 9 n = 0, where n is the normal to the exit plane. This con- 
dition is generally valid if the extent of the computational domain in the 

primary flow direction is sufficient to ensure fully-developed flow conditions 
for internal flows or self-similarity for jet flows at the exit plane. Particle 
trajectory and temperature history computations are discontinued for a compu- 
tational particle once the particle exits from the computational domain through 

the exit plane or any other open boundary. 

Solid Boundary: 

The conventional wall functions approach is used to impose wall boundary 
conditions on the velocity and temperature as well as the turbulence kinetic 
energy and its dissipation rate. The presence of particles in a fluid flow has 
been experimentally observed to influence the boundary layer [Kramer and Depew, 
1972] and, as a consequence, the nature of the wall function which is normally 
used to connect the actual value of a given variable at the wall to the value 
at the wall-adjacent grid node. During their trajectories, particles that reach 
the wall either adhere to it as observed in particle erosion problems [Dosanjh 
and Humphrey, 1984], or collide with the wall and get "reflected” back into the 
flow domain, usually with an accompanying loss of energy and momentum to the 
wall. In addition, the high level of shear in the wall vicinity coupled with a 
particle velocity slip introduces an additional transverse force on the par- 
ticle which further modifies its subsequent trajectory and behavior in the 
near-wall region. These effects have not been included in the present study, in 
which perfectly reflecting boundary conditions have been adopted for the parti- 
cle-wall interaction, but will be the subject of a future study. 

Other generic boundary condition types include the symmetry axis, for 
which /5n = 0, where in this case, n is the normal to the symmetry axis, 
and the open boundary condition which has been used by Leschziner and Rodi 
[1984], Dosanjh and Humphrey [1984], Amano and Brandt [1984] and Chen and Ade- 
niji-Fashola [1987] for modeling parabolic flows of free jets and wall jets using 
elliptic formulations. These are described in greater detail by Adeniji-Fashola et al. 
[1988]. 
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6. EXAMPLE FLOWS 


In order to illustrate the multiple-realization particle trajectory mode- 
ling scheme for confined turbulent fluid-particle flows described above, two 
example flow problems - vertical pipe upflow and horizontal coaxial jet flow in 
a recirculation chamber with a particle-laden central jet and a clean annular 
jet are examined. 1500 computational particles were found to be adequate in 
each example for obtaining statistically stationary solutions. Typically, 
global under-relaxation values of 0.50 were found adequate to ensure the sta- 
bility of global iterations of which between five and seven were required to 
obtain globally converged solutions. The results obtained for the numerical 
simulation of these flows will now be discussed. 


VERTICAL PIPE UPFLOW 

The experimental data which served as the basis for this example numeri- 
cal simulation are those of Tsuji et al. [1984] for the upflow of a particle- 
laden stream in a straight vertical pipe. The experimental flow within the test 
section is considered to be fully-developed after going through a riser that is 

167.5 diameters long. 

A 50 X 23 uniform grid distribution was used to discretize the computational 
domain which had an axial extent of 60 pipe diameters. Figure 3 shows both the 

experimental data and the numerical predictions of the radial profile of the slip in 
the axial velocity between the air and the particulate phase. The mean particle size 
and loading ratio are 200 pm and 1.0 respectively. The air velocity is slightly 
overpredicted in the 0.2R - 0.8R range where R is the pipe radius. However, the 
prediction accuracy is considered to be good for such a complex system. The radial 
profile of the axial velocity of the solid phase is particularly well predicted. The 
location of the cross-over in sign of the slip between the phases is predicted to be 
closer to the wall, less than 0.1R from the wall, than the 0.2R from the wall that 
was experimentally observed. 

A similar picture obtained for the higher loading ratio of 2.1 is pre- 
sented in Figure 4. The level of accuracy of the predictions is similar to that 

of the 1.0 loading ratio case. However, it is the air velocity profile that is 

better predicted in this case. The solid phase axial velocity is considerably 
underpredicted in the inner 60 percent of the wall region. 

As pointed out earlier, the particulate phase has the effect of modulating 
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the level of the turbulence intensity. For smaller particle sizes this results 

in a decrease in the kinetic energy of turbulence. The experimentally observed 
and numerically predicted turbulence modulation effect for a loading ratio of 
3.2 are illustrated in Figure 5. The solid line in the figure shows the pre- 

dicted radial profile of the turbulence intensity for the corresponding "clean" 
flow. The predicted level of turbulence intensity is considerably higher than 

the level observed from experiment. Also, while a greater modulation effect was 

observed closer to the wall region, the predictions show a reversal in which 
the greater level of modulation is located closer to the pipe centerline. The 
imposed wall boundary conditions and wall functions in the numerical scheme are 
probably responsible for the suppression of the modulation effect in the 

near-wall region. 

The development in the axial direction of the streamwise velocity of the 
particulate phase for an inlet velocity slip ratio of 0.10 is shown in the 

contour plot of Figure 6a and a corresponding surface plot in Figure 6b. The 
ability of the particle trajectory scheme to effectively handle extreme levels 

of velocity slip was tested by imposing an axial slip velocity of 0.10 at the 
pipe inlet plane. The figures indicate that a fully developed state was 

attained in the 60D extent of the computational domain. 


HORIZONTAL COAXIAL JET FLOW IN RECIRCULATION CHAMBER 

The experimental data of Leavitt [1980] serve as the basis for the numer- 
ical simulation of this example. The actual geometry studied is illustrated in 
the schematic of Figure 7. The primary jet air velocity at inlet is 33 m/s 
while the corresponding secondary jet air velocity is 42 m/s. Coal particles of 
a mass mean diameter of 43j*m were used to uniformly seed the primary jet and 
the particle loading ratio is 1.50. The estimated turbulence intensity levels 

at the inlet are 15 and 18% for the primary and secondary jets respectively. 

The primary and secondary jet diameters at inlet are 0.0255m and 0.127m respec- 

tively while the chamber diameter is 0.206m. The axial extent of the recircula- 
tion chamber is 0.926m (36.3 primary jet diameters or 4.5 chamber diameters). 

A 41 X 41 non-uniform staggered grid distribution, shown in Figure 8, is 
used for the numerical study and the computational domain extended to 20D where 
D is the chamber diameter. The numerical prediction of the evolution of the 
axial velocity is shown in Figure 9. The corner recirculation zone is seen to 
extend to about 1.79D. No particles are predicted as reaching this recircula- 
tion zone and this is believed to be due to the high chamber-to-primary jet 
diameter ratio of 8.08 and the positive slope of the shear in the mixing layer 
between the primary and the secondary jets which will result in a slip-shear 

transverse force directed towards the centerline. Another interesting observa- 
tion is that the particle axial velocity starts lo lead that of the fluid from 
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about the 0.80D axial location and this continues to about the 7.43D axial 

location downstream of which all axial velocity slip disappears. Particles are 

seen to have dispersed to the outer extremities of the recirculation chamber by 

the time the 12.0D axial location is reached. However, it should be remembered 

that this is only a hypothetical situation since the actual experimental inves- 
tigation was limited to an axial extent of only 4.5D. 

Figure 10 shows the axial evolution of the turbulence intensity. It is 
observed that up to about the 3.0D axial location, the turbulence intensity in 

the presence of particles (shown dotted) falls below that of the clean flow in 
the primary jet portion of the flow but is actually higher for the rest of the 

chamber in the radial direction. However, beyond the 3.0D axial location, the 

clean flow turbulence intensity uniformly lags the two-phase intensity at all 

radial locations for any given axial location. The kinetic energy of turbulence 
is essentially fully developed at the 5.15D axial location and only a radially 
uniform decrease in magnitude is observed for the rest of the flow in the axial 
direction. This is in contrast to the radial profile of the axial velocity 
which does not become fully developed for both phases until the 12.0 to 15.0 
diameter range is reached. 

The contour and surface plots of the particle axial velocity are shown in 
Figures 11a and lib. These have been normalized with respect to the secondary 
jet gas velocity at inlet. Since, in contrast to the two-fluid scheme, non-zero 
values of the particle velocity are not returned for computational cells not 
visited by any particle during the trajectory calculations, the zero-velocity 
surface in the plots of Figure lib also indicate the particle-deficient 
regions. 

The comparison of the limited experimental data available from Leavitt 
[1980] is currently being undertaken. 


7. CONCLUDING REMARKS AND RECOMMENDATIONS FOR FUTURE WORK 


The numerical modeling of confined turbulent fluid-particle flows using 
the multiple-realization particle trajectory scheme has been presented. The 
performance of the numerical modeling scheme has been tested using data for the 
upward flow of a fluid-particle stream in a straight vertical pipe and for the 
horizontal coaxial jet flow in a large recirculation chamber for which the cen- 
tral jet is particle-laden. 

The multiple-realization particle trajectory turbulent flow modeling 
scheme ... 
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is believed to be a more elegant and accurate approach to the 
extension of single-particle hydrodynamics to dilute multi-particle 
systems; 

is better able to incorporate additional force terms such as lift, 
virtual mass and Bassett history terms in the particle equation of 
motion as appropriate; 

needs further investigation in order to improve its computational 
efficiency and so reduce its huge CPU time requirements; 

needs to have the particle-turbulence and particle-wall interactions 
further investigated to improve prediction accuracy. 
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